Level-statistics of Disordered Systems: a single parametric formulation 



Pragya Shukla* 

(i) Max-Planck Institute for Physics of Complex Systems, Nothnitzer Str.38, Dresden-Oil 87, Germany, 
(ii) Department of Physics, Indian Institute of Technology, Kharagpur-721302, India. 

(February 1, 2008) 

We find that the statistics of levels undergoing metal-insulator transition in systems with multi- 
parametric Gaussian disorders behaves in a way similar to that of the single parametric Brownian 
ensembles The latter appear during a Poisson — > Wigner-Dyson transition, driven by a ran- 
dom perturbation. The analogy provides the analytical evidence for the single parameter scaling 
behaviour in disordered systems as well as a tool to obtain the level-correlations at the critical point 
for a wide range of disorders. 
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The nature of the eigenfunctions can significantly af- 
fect the statistical behaviour of energy levels of a dis- 
ordered system and thereby related physical properties 
e.g transport. The presence of disorder may cause lo- 
calized waves in the system, implying lack of interaction 
between certain parts. This is reflected in the structure 
of the Hamiltonian matrix which is sparse or banded in 
the site representation depending on the dimensionality 
of the system. The variation of the disorder-strength can 
lead to a metal-insulator transition (MIT), with eigen- 
fuctions changing from a fully extended state (metal) to 
a strongly localized one (insulator) with partial localiza- 
tion in the critical region. The associated Hamiltonian 
also undergoes a transition from a full matrix to a sparse 
or banded form and finally to a diagonal matrix. The 
statistical studies of levels for various types of disorders 
require, therefore, analysis of different ensembles. Here 
the nature of the localization and its strength is reflected 
in the measure and the sparcity of the ensemble, respec- 
tively. Our objective in this paper is to obtain a math- 
ematical formulation for the level-correlations, common 
to a large class of disorders (Gaussian type); the informa- 
tion about the nature of disorder enters in the formula- 
tion through a parameter, basically a function of various 
parameters influencing the localization. 

Recently it was shown that the eigenvalue distributions 
of various ensembles, with a Gaussian measure, appear 
as non-equilibrium stages of a Brownian type diffusion 
process ||. Here the eigenvalues evolve with respect to 
a parameter related to the complexity of the system rep- 
resented by the ensemble. The solution of the diffusion 
equation for a given value of the parameter gives, there- 
fore, the distribution of the eigenvalues, and thereby their 
correlations, for the corresponding system. The present 
study uses the technique to analyze the spectral proper- 
ties of the levels undergoing MIT in various disordered 
systems. 

The Anderson model for a disordered system is 



described by a d-dimensional disordered lattice, of 
size L, with a Hamiltonian H = ^ n e n a^a n — 
Hn^ra & nm( a n a m + a„a+) in tight-binding approxima- 
tion. The site energies e n , measured in units of the over- 
lap integral between adjacent sites, correspond to the 
random potential. The hopping is assumed to connect 
only the z nearest-neighbors (referred by m) of each site. 
In the site representation, H turns out to be a sparse 
matrix of size N — L d with diagonal matrix elements 
Hkk — £fc and off-diagonals Hki describing the hopping. 
The level-statistics can therefore be studied by analyzing 
the properties of an ensemble of (i) sparse real symmetric 
matrices, in presence of a time-reversal symmetry and (ii) 
sparse complex Hcrmitian matrices in absence of a time- 
reversal. 

The MIT can be brought about by a competitive varia- 
tion of the disorder and hopping rate which can be mim- 
icked by a change of the distribution parameters of the 
ensemble measure p{H). In this paper, we consider the 
case in which the site-energies e, are independent Gaus- 
sian distributions. The hopping can be chosen to be non- 
random or random (Gaussian). The p(H), for any inter- 
mediate state of MIT can therefore be described by 



p(H, y, b) = Cexp[- ^ £)(l/2/i w . a )(fl*,. a ~ bkl ^ W 

8=1 k<l 

with subscript "s" of a variable referring to its compo- 
nents, (3 as their total number {(3—1 for real variable, 
(3 = 2 for the complex one), C as the normalization con- 
stant, h as the set of the variances hu-,s =< H^. s > 
and b as the set of all bki- s - As obvious, in the limit 
hki-i,hki-,2 0, cc l-(l) corresponds to the non-random 
nature of H k i {H k i = b k ui + ibu^)- Here the strongly 
insulated state corresponds to all hki — ► (k ^ I) and 
bki —> (all k, I), implying no overlap between site states. 
The metal, with almost similar overlap between various 
site states, can be modeled by eq.(l) with all hki — * 7" 1 , 
that is, a Wigner Dyson (WD) ensemble pi. 
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The eigenvalue distribution for a metal is given by 
the WD distribution, and, for an insulator by a Pois- 
son distribution The distribution for various tran- 
sition stages can be obtained by integrating p over the 
associated eigenvector space. Let P(p, h, b) be the prob- 
ability of finding eigenvalues A^ of H between fii and 
fii + d/ii at a given h and b, it can be expressed as 
P(/i, h,b) = J rii=i ^(m* — ^i)p{H, h, b)dH. As discussed 
in ref. Q , it is possible to define a " complexity" parame- 
ter Y, a function of various distribution parameters hu- s 
and b k i- s 0, 
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(here the Yl implies a product over non-zero bki-, s and 
Xki-s = 1 - ^guhkUs only, gui = 2 - fe, C as a constant 
determined by initial ensemble and M as the number of 
all non-zero parameters xm- s and bki-, s ) such that the evo- 
lution of P with respect to Y results in the diffusion of 
eigenvalues with a finite drift due to their mutual repul- 
sion, 
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(f3 — 1,2 for Hamiltonians with and without time- 
reversal, respectively). Here 7 is an arbitray param- 
eter M. The evolution reaches a steady state when 
dP/dY — ► 0; in this limit is given by a WD dis- 

tribution, P(ji) = H t<j \m - fijfe-l ^ fc 

The spectral correlations for MIT can now be obtained 
by using following analogy. The eq.(3) is same as the 
equation governing an evolution of the distribution P of 
the eigenvalues, of a Hamiltonian H = ^/J(Hq + XV) 



with A 2 



3 2 7 (F-Yo) 



1 and / = e- 2 ^ Y ~ Y °\ from 



initial state Hq to WD ensemble. The transition, re- 
ferred as WDT later on, is caused by a random per- 
turbation V, taken from a WD ensemble (of variance 
7), and of strength A Jl],||]; the elements H^i are thus 
Gaussian distributed with a variance h^i = (1 — f)/2j, 
hkk = 1/27 and same mean for all of them. The tran- 
sition (WDT) to equilibrium, with Y — Yq as the evo- 
lution parameter, is rapid, discontinuous for infinite di- 
mensions of matrices But for small- F and large N, 
a smooth crossover can be seen in terms of a rescaled 
parameter A = (Y — Y )/A 2 , with A,, = AN/tj as 
the mean-level spacing in the correlated region of the 
spectrum at Y — Yq, r\ as the correlated "volume" and 
A = A(/x, Y) as the mean level spacing of the whole 
spectrum. The intermediate states, corresponding to fi- 
nite, non-zero A-values, during the crossover are known 
as Brownian ensembles (BE) ]l],|J or Rosenzweig-Porter 
ensembles (RPE) The similar evolution equa- 

tions of P for MIT and WDT, imply a similarity in their 



eigenvalues distributions for all F-values and thereby 
correlations for all A-values, under similar initial con- 
ditions (that is, P(/j,,Yq) same for both the cases al- 
though p(H,Yo) may be different). In finite disordered 
systems, therefore, a continuous family of intermediate 
statistics exists between metal and insulator which can 
be described by the BE for the corresponding A-value, 
occuring during Poisson — > WD ensemble type WDT. 

The behaviour of the mean level spacing A and the cor- 
relation volume X] divide the level-statistics for the large 
BE (N — > 00) into three regions 0: 

(i) N 2 (Y - Y ) -> 0: A cx TV" 1 and r) < N, giving 
A — » and Poisson statistics, 

(ii) N 2 (Y - Y Q ) -> 00: A cx TV -1 / 2 and rj > N giving 
A — > 00 and WD statistics. 

(iii) A 2 (y — Yq) = l/2c, with c as an arbitrary con- 
stant with respect to TV: although A still behaves as 
o(l /N), but now 77 ss TV thus giving A = l/2c7r (referred 
as A. be)- 

The A = A* BE therefore corresponds to a third statis- 
tics, intermediate between Poisson and WD ensemble 
and is known as the critical Brownian ensemble (CBE). 
This being the case for arbitrary values of c, an infi- 
nite family of CBE occur during WDT. The presence 
of such a family can be seen from any of the fluctuation 
measures for WDT. One traditionaly used meausure in 
this regard is relative behaviour of the tail of nearest- 
neighbour spacing distribution P(s), defined as a(S, A) = 
J Q S {P(s) - P w (s))ds/ Jq(P p (s) - P w (s))ds with 6 as any 
one of the crossing points of P w (s) and P p (s) (here sub- 
script w and p refer to the WD case and Poisson case 
respectively). In the limit TV — » 00, a = and 1 for 
WDE and Poisson limit respectively. The FIG.l shows 
the numerically obtained behaviour of a with respect to 
|c — c*\ (for 5 w 2.02), with c* corresponding to one of 
the CBEs. The convergence of all the points on two 
branches for different A-values confirms the existence 
of a CBE at c* with a critical exponent v — > 00 (as 
a(c) w a{c*) + constant.\c~ c*|A 1 / l/ near c*). The TV- 
independence of ot(c) also indicates the possibility of in- 
finitely many c* (and its arbitrarity) and therefore an 
infinite family of CBEs. 

For disordered systems with infinite system size L, the 
states are critical near a particular disorder W* (or en- 
ergy), termed as critical point, with a correlation length 
G(W) cx \W — W*\~", v as the critical exponent. At 
the critical point, £ sa L and the critical value of the 
parameter A (referred as A* AH ) depends on the dimen- 
sionality; A* AH = (Y - Y )/A 2 , with both Y and A, ; 
dimensionality-dependent. A knowledge of A* AH can then 
be used to map the critical level statistics at MIT for 
various dimensions d > 2 — ► 00 to the infinite family of 
critical Brownian ensembles (CBE). 

The identification of the appropriate CBE correspond- 
ing to the critical Anderson Hamiltonian (CAH) rc- 
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quires h* BE = A* AH , and, thus a prior knowledge of A- 
value associated with CAH. We consider here one ex- 
ample (this case is also used in our numerical anal- 
ysis). Consider an Anderson system with the Gaus- 
sian disorder, same for each site, and random or non- 
random hopping between nearest neighbours. The cor- 
responding ensemble measure can be described by eq.(l) 
with h kk = W 2 /2, b kk = and h kl = W? /2, b u = 
t for {fc, 1} pairs representing hopping, h k i — > and 
b k i — > for all {k, 1} values corresponding to diss- 
connected sites. This gives, by using eq.(2), Y — 
-{N/2M~f)lTi[\l - jW 2 \\l - 2~ f W?\ f3z / 2 \t + Stol 13 *] + C. 
Here M = (3N(N + z5 t0 + 2 - j3)/2 with zN as 
the number of the connected sites (nearest-neighbours) 
which depends on the topology and the dimensional- 
ity d of the system. Analogous to WDT, a rescal- 
ing of Y — Yq would give the parameter for smooth 
transition: A = 2(a — ao) /3 _1 (C / ' L) d with a(a,t) = 
ln[|l-7^ 2 ||l-2 7 M^ 1 2 |' 32 / 2 |t + (5 t o| /32 ], and a as the 
value of a at Yq. Here A behaves as 0(N~ 1 ^ 2 ) (as 
A^ 2 (Y" ~ Yq) — > oo) and -q — Q d with ( as the localiza- 
tion length which depends on disorder W and can be 
determined by using wavefunction statistics (e.g. in- 
verse participation ratio). As disorder decreases, the 
C(W) increases resulting in a smooth increase in A in 
finite systems, and, thereby the statistics intermediate 
between Poisson and WDE. However for infinite system 
sizes, A — * and the statistics is Poisson as long as 
£ < L. But at a certain disorder strength, ( ~ L, which 
makes A size-independent thus giving its critical value 
A* AH w 2(a — ao)/? -1 . The level statistics of the sys- 
tem at this disorder strength is given by the CBE with 
the same A-value and is critical due to (i) being different 
from both Poisson and WD behaviour, (ii) its invariance 
as N — > oo. Now for a further decrease in disorder, £ > L, 
as a result A — * oo, and the system has a WD statistics. 
However, for an infinitely long one dimensional lattice 
with z nearest neghbours, the statistics always remains 
Poission irrespective of disorder strength. This is because 
£ w z 2 giving A z 2 /L — > for L — > oo unless there is a 
very long range connectivity in the lattice (i.e z ~ VL). 
However a crossover from Poisson to WD ensemble can 
be seen for finite L by varying the ratio z 2 /L. 

In general, A will be a function of coordination num- 
ber, disorder strength, hopping rate and dimensionality 
of the lattice as well as the level-density of the system. 
The boundary conditions/ topologies, leading to different 
sparcity and coordination numbers, may therefore result 
in different critical level statistics even if degree of dis- 
order, hopping rate and dimensionality is same; this is 
in agreement with numerical observations [fill and ana- 
lytical study for 2D systems |l2]| . Similarly different di- 
mensions, can lead to diffenent local level-densities and 
thereby different critical points A*. 

Many results for the spectral fluctuations of the WDT 



with Poisson ensemble as an initial state are already 
known ^ and can directly be used for the correspond- 
ing measures for the MIT in different disordered sys- 
tems. For example, consider the 2-level density correla- 
tor R 2 (r;A); R 2 (r;A) =< v^, A)v([j, 2 , A) > / < v > 2 = 
/ P(n, A)d M3 ..dMJV. Here vQt, A) = A^ 1 J2 t S(ji- 
l«i) is the density of eigenvalues and < .. > implying the 
ensemble average. The R 2 for the Anderson transition 
in presence of a magnetic field, can therefore be given 
by the R 2 for WDT between Poisson — > GUE ensemble 
|,|; R 2 = l-Y 2 where 

Y 2 (r; A) = -— / duF e- 2Au2 ~ 47rAu (4) 
r Jo 

with F = sin(ur)/i — cos(ur)/ 2 , f\ — (2/z)[Ii(z) — 
y g7jjja0 0]. h = (l/u)[h(z) - ^/2^I 3 (z)}, z = 
V327rA 2 u 3 and /„ as the n th Bessel function. Here 
R 2 (r, oo) = 1 — (sin 2 (7rr)/7r 2 r 2 ) and R 2 (r,0) = 1 corre- 
sponding to metal and insulator regime respectively. A 
substitution of A = A* in eq. (4) will thus give the R 2 for 
the CAH. Similarly the nearest-neighbour spacing distri- 
bution P(s) for the MIT can be given by using the one for 
the BE during Poisson — > GUE transition P(s; A) oc 
(2^A)- 1 /2 se - s2 / 8A r °° dx x- 1 e- x - x ^ SA sinh(xs/4A). 
Further the good agreement of the numerically obtained 
P(s) for the CAH without time-reversal symmetry with 
that of a CBE (with same A-value) reconfirms our an- 
alytical result; the FIG. 2 (a, b) show the behaviour of 
the CAH in presence of a random hopping (CAH-I with 
W = 8.15, W x = 1, t = 0, and, thereby A = 10.6) 
and a non- Random hopping (CAH-II with W — 21.3, 
W\ = 0, t = 1, thus A w 0.53), respectively along 
with their CBE analogs (c determined by the relation 
A* AH = A* BE = (2c7r) _1 ); see § for the details on the 
numerics. 

An important characterstic of critical level statistics is 
the level compressibility Xj x(A) ~ 1 — /_ l2(^;A)dr; 
X = 0, 1 in the metallic and the insulator phase, re- 
spectively, and takes an intermediate value at the hybrid 
phase near the critical point. The eq.(4) can be used to 
obtain x(A) = 1 - 47rA du/i(z)exp[-2Aw 2 - 4nAu]. 
The critical region, with its finite A value (= A*), will 
thus have a level compressibility different from both 
metal and insulator regimes; a < % < 1 value is sup- 
posed to be an indicator of the multifractal nature of the 
eigenvectors. 

The compressibility of the spectrum can also be seen 
from the " number variance" Ti 2 which describes the vari- 
ance in the number of levels in an interval of r mean level 
spacings. The Tj 2 (t) : basically a measure of the "spec- 
tral rigidity" is related to the compressibility: lim r — > 
oo S 2 ~ x r - The FIG. 3 shows the numerical behaviour 
for the E 2 (r;A) for the CAH-I and CAH-II along with 
their CBE analogs; the good agreement in each case re- 
verifies our claim about the similarity between the MIT 
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and WDT. Furthermore FIG. 3. also shows the fractional 
behaviour of x f° r CBE, which is in agreement with the 
X-form given in the preceeding paragraph. 

The statistical measures for the Anderson tran- 
sition in presence of a time-reversal symmetry can 
similarly be obtained by using their equivalence to 
WDT from Poisson — > GOE ensemble. However due 
to the technical difficulties @]), only some approxi- 
mate results are known for the latter case, for ex- 
ample, the i?2 for small- r can be given as Ri(r, A) w 
(2r - 1) 1/2 J1/3 ((2r - l) 3 / 2 /3A) e r / 2A with J(z) as the 
Bessel function (obtained by solving eq.(17) of ||). Sim- 
ilarly for large-r behaviour, i? 2 can be shown to satisfy 
the relation (see eq.(21) of R 2 (r,A) = R 2 (r, oo) + 
2/3 A dsi? 2 (r - s; 0) - i? 2 (r - s; oo)/(.s 2 + 4tt 2 /? 2 A 2 ). 
By taking 13 — 1 and i? 2 (r, oo) = 1 — sin 2 (7rr)/7r 2 r 2 — 
( J* 00 dxsin7rx/7ra;) (^sinTrr/Trr 1 ) (GOE limit), the 
R 2 (r, A) for the transition in presence of a TRS can 
be obtained; i? 2 (r, A) w # 2 (r, oo) + 4A/(r 2 + 4tt 2 /3 2 A 2 ) . 
The lack of the knowledge of i? 2 (r, A) for entire energy- 
range handicaps us in providing an exact form of 
However the P(s) for this case can be given by using the 
one for a BE (RPE) during Poisson -^GOE transition 



@: P(s,A) = (tt/SA) 1 / 2 ^ 
as the Bessel function. 

In the end, a comparison of our results with some of 
the past studies is crucial. The presence of a fractional 
compressibility in an ensemble essentially same as CBE 
and its possibility as a model for CAH was also sug- 



gested in|9] which was later on contradicted in ]lC| |. The 
claim in llC ] is disproved by our numerical studies of both 
the transitions, showing similar behaviour for compress- 
ibility, besides other fluctuation measures at the critical 
point. The numerical work therefore reaffirms our claim 



based on the exact analytical work: (i) the level-statistics 
for MIT in disordered systems with a Gaussian disorder 
can be described by the same for WDT with a Poisson 
initial condition, (ii) the level-statistics in the disordered 
systems is indeed governed by a single scaling parameter. 
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I. FIGURE CAPTION 



CAPTIONS 

Fig. 1. a vs |c- c*\ for WDT 

Fig. 2. P(S) vs S for (A) CAH-I and the corresponding CBE (c = 0.015), (b) CAH-II and the CBE with c = 0.3. 
Fig. 3. E 2 vs r for (i) CAH-I and CBE with c = 0.3, (ii) CAH-II and the CBE with c = 0.015. 
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